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We propose a model for the dynamics of a social system, which includes diffusive effects and a 
biased rule for spin-flips, reproducing the effect of strategic choices. This model is able to mimic 
some phenomena taking place during marketing or political campaigns. Using a cost function based 
on the Ising model deflned on the typical quenched interaction environments for social systems 
(Erdos-Renyi graph, small- world and scale-free networks), we flnd, by numerical simulations, that a 
stable stationary state is reached, and we compare the final state to the one obtained with standard 
dynamics, by means of total magnetization and magnetic susceptibility. Our results show that the 
diffusive strategic dynamics features a critical interaction parameter strictly lower than the standard 
one. We discuss the relevance of our flndings in social systems. 



I. INTRODUCTION 

In the past few years the application of statistical me- 
chanics to social phenomena gave rise to interesting mod- 
els, which were able to capture some general mechanisms 
in opinion forming. In these models, the relations be- 
tween people in a group are represented by a network 
with a given topology, where sites are people, links model 
interactions and the opinion of a single agent is typically 
represented by a discrete variable on the corresponding 
site. 

One of the most important aspects in opinion forming 
within a community is the dynamics through which de- 
cisions of single agents take place, possibly leading the 
system to a stationary final state. In this framework 
different approaches have been considered, in order to 
investigate the configuration reached by the system at 
different time scales. A common approach, as in voter 
model and Axelrod dynamics (see e.g. [H-Q), is to intro- 
duce directly a dynamical rule according to which dis- 
crete variables, representing the choice of a single agent, 
evolve. Then, the average value of these discrete variables 
and its dependence on the dynamics itself as well as on 
the initial conditions is investigated. However, in many 
situations, ranging from polls to marketing analysis, rel- 
evant global parameters describing the behavior of large 
but finite subsamples of the populations are not rapidly 
changing on the time scale considered. If the results of 
the experimental measurements are stable, then it can 
be meaningful to analyze the social system by looking to 
its equilibrium or stationary behavior. This approach, 
directly linked to statistical mechanics and used in so- 
cial sciences, is based on the introduction of a dynamics 
derived from a cost function H 043 j which depends on 



the configuration of the system and on a set of param- 
eters, measuring the interaction between people in the 
community. The cost function theory has also shown 
its potentials in micro-economy especially thanks to the 
work of Mc Fadden [8]. There, quantitatively precise 
predictions were made on social behavior by means of 
discrete choice theory, a simple probabilistic approach 
based on the individual independent choice assumption, 
which corresponds to non-interacting agents. Among the 
main purpose of the present paper there is the attempt to 
extend to interacting systems the approach pioneered in 
fS] and successfully applied to social systems by Durlauf 
4] within the framework of equilibrium statistical me- 
chanics of deterministic systems. Our work here ex- 
tends those previous findings in several directions: first 
it assumes that there is interaction (including the ran- 
dom case) among the agents and investigate the possi- 
ble steady states by means of natural dynamical systems 
usually encountered in social behavior. The dynamics 
is implemented with a cost function in which the pa- 
rameters measuring the interactions between people can 
be assumed to be frozen. In real world phenomena the 
dynamics has two built-in time scales: one for the opin- 
ion flip, the other for the interaction change. In most 
of the social and economical examples the interactions 
among individuals change at a rate enormously slower 
than the opinions of each of them. In modeling such sit- 
uations the interactions variables will be considered as 
frozen during the whole run of the dynamics and the 
physical quantities will be averages of stationary states. 
The procedure is known as "quenched" stationary state 
and parallels the standard quenched equilibrium measure 
of disordered systems (like spin glasses) in which the free 
energy is computed averaging in the disorder after the 
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logarithm. 

Once the cost function and the related evolution dy- 
namics are identified, the system can be studied numeri- 
cally, for example via Monte Carlo simulations. In partic- 
ular, if one expects that the final state reached by the sys- 
tem will obey an equilibrium Boltzmann distribution, the 
dynamics is usually implemented to lead the system to 
that state, imposing a detailed balance condition. How- 
ever, this is quite unrealistic in real social communities. 

As a matter of fact, dynamical evolution within social 
networks often exhibits a strategic character: each agent 
often chooses between his neighbors and decides which 
link to activate at a given time, according to his own 
advantage, based on the observation of what happens to 
his neighbourhood. 

A common situation is the following: people try to con- 
vince their neighbours in order to share the same opin- 
ion or trait because of ideological reasons or, also, be- 
cause this translates in some economical advantage. To 
fix ideas let us suppose that the opinion or the trait con- 
sidered can be described by a binary variable. For exam- 
ple, during an electoral campaign or before a referendum, 
people try to convince their acquaintances to support a 
given candidate or a given position. In a different con- 
text we can think about a community of people where 
each agent has made a subscription to a phone company. 
Let us suppose that only two different companies, A and 
B, exist so that we can distinguish between A- users and 
B-users. Now, fares for phone-calls are different accord- 
ing to whether the call occurs between customers of the 
same company or between customers of different com- 
panies, being higher in the latter case. As a result, for 
an A-user (B-user) the optimal situation is when all his 
acquaintances are also A-users (B-users) as he can then 
enjoy low fares. Hence, each agent would like to induce 
his neighbours to adopt the same company. 

In both situations cited above, strategies are possible if 
an agent knows the neighbourhood of his acquaintances. 
As an example, if I want to increase the number of, say, 
A-users, among my acquaintances, I can either pick up 
a friend of mine randomly among those who are B-users 
and try to convince him to become an A-user, or select 
among my friends the B-user whose acquaintances are 
mostly A-users. The latter strategy is of course expected 
to be more effective. 

Another important aspect in dynamics concerns the 
rule according to which agents are selected and given the 
possibility to change opinion. Of course, a deterministic 
updating, though computationally efficient, is unrealis- 
tic: the dynamics must contain a stochastic character, 
reproducing the randomness typical of social interactions 
[ol-fl^. More likely, we can assume that the opinion up- 
dating results from phone calls, mail exchanges and other 
contacts among agents and this exhibits a diffusive fea- 
ture. Hence, the newest updated agent will choose among 
his neighbours the next agent to be updated. 

In this paper, we intend to study the equilibrium 
reached by the system, where an explicit cost function 



has been introduced, endowed with a dynamics which 
takes into account these two aspects: strategy and diffu- 
sion. We first focus on the dynamical process, modelling 
the quenched social network by an Erdos-Renyi random 
graph. This graph provides a stochastic network, able to 
capture some aspects of a real community, and allowing 
for some exact calculations (T3l - [l5| . We then extend our 
results to scale free and small- world networks, which are 
known to reproduce the typical topological features of 
real social networks [iMI. We adopt as a cost function 
the ferromagnetic Ising Hamiltonian which, being one of 
the simplest model mimicking interactions amongst the 
agents of social systems, allows us to focus on the dynam- 
ical process. In particular, the interaction parameter J 
here represents the "imitation strength" and it measures 
how important it is for two nearest-neighbours to agree. 
For instance, in the phone companies example, a large 
value of J corresponds to a situation where a phone call 
between A-A or B-B users is much cheaper than a phone 
call between A-B users. We focus here on ferromagnetic 
interactions as these are considered to be the predomi- 
nant feature of social interactions in several contexts [20| . 
Our dynamics can however be extended to the case of 
more complex interaction patterns. 

We find that after a suitable relaxation time the values 
of the global observables of the system display time aver- 
ages independent of the initial conditions and whose fluc- 
tuations decrease with the system size, indicating that a 
stationary situation is reached. We also recover the phase 
diagram expected for the Ising model on the Erdos-Renyi 
random graph, small-world and scale-free networks, and 
this implies that in all these cases there exists an inter- 
action parameter Jc such that if J < Jc the number of 
A-users equals, on the average, the number of B-users, 
while if J > Jc a prevailing group emerges. However, 
we evidence a remarkable difference: With respect to the 
case of a non-strategic dynamics, the critical region is 
always shifted to a lower value of the interaction param- 
eter J. In the example of the competition between the 
two phone companies, this means that, once the price 
policy has been set by the companies, i.e. once a given 
J > Jc has been fixed by the price differences, if the 
equilibrium is reached by a strategic dynamics, the ex- 
tent of the prevailing community is larger. The existence 
and the position of a transition point represent a key in- 
formation in a social system, as they signal an unstable 
situation, which is to be favored or avoided, depending 
on the meaning of the global parameter. Moreover, it is 
important to understand what influences the position of 
the critical region in order to control and tune it, if pos- 
sible. Our result suggests that the onset of the critical 
regime is indeed influenced by the dynamics. 

The paper is organized as follows: in Sec.|ll]we present 
our model and we describe the topologies it is embedded 
in. Then, Sec. HIT] is devoted to the description of the 
strategic diffusive dynamics introduced and in Sec. II VI we 
show our results. Finally, Sec. |V] is left for conclusions 
and outlook. 
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II. MODEL AND NOTATIONS 

A social network is meant as a (typically large) set 
of people or groups of people, also called "agents" , with 
some pattern of interactions between them. This can be 
efhciently envisaged by means of a graph whose nodes 
represent agents and links between two of them repre- 
sent the existence of a relationship (which could be ac- 
quaintanceship, friendship, etc.). Therefore, each agent 
i is connected with a set of "nearest-neighbours" , whose 
number is referred to as the "degree" of the node i. 

Now, several kinds of graph have been proposed in the 
past as models able to mimic the features displayed by 
a real population, and they are all built starting from 
three main topologies: random graphs, small- world and 
scale-free networks. The random graph introduced by 
Erdos and Renyi (ER) [2lj is one of the most studied 
since it combines a stochastic character with an easy def- 
inition which allows to calculate exactly many interesting 
quantities [l5|. However, real social networks have been 
shown to feature some peculiar topological and metric 
properties, which are not all included in the ER random 
graph. Two typical features are: power-law degree distri- 
bution (scale- free topology), shortness of geodetical path 
(small- world phenomenon) . 

The ER random graph can be defined as follows: given 
a number TV of nodes, we introduce connections between 
them in such a way that each pair of vertices i,j has 
a connecting link with independent probability p. As a 
consequence, the probability Pa that a node in a random 
graph has degree exactly equal to a is given by the bi- 
nomial distribution 
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The small-world network can 
be built starting from a ring with N nodes, and "su- 
perposing" to the ring a random graph with given mean 
connectivity a, by adding a link to two points i and j on 
the ring with probability pa = a/N (but it is possible to 
obtain small- world topologies with different rules [2^). 
To account for the large variability of the degree in real 
social networks, one of the most studied topologies is the 
scale-free graph, where the probability Pa that a node 
has degree equal to a is given by a power law: 



Pa 



a 



(1) 



with 1 < 7 < 3. 

We now outline the general framework for modeling in- 
teractions among agents. First of all, we associate to each 
agent i a binary variable Si ~ ±1, representing the two 
possible forms of the considered opinion or trait. For ex- 
ample. Si = might indicate that the i-th agent does 
support the current government or is an A-user, while 
Si = —1 that he does not support the current govern- 
ment or that he is a B-user. The whole community, de- 
scribed by the set s = {si, S2, sat}, will therefore be 
characterized by the mean value 



which can be measured by, say, a referendum vote or a 
survey. 

We assume that agents do not possess any a priori bias 
towards +1 or —1 state, but they move towards a given 
trait as a result of the interaction with their nearest- 
neighbours. More precisely, we introduce a cost function 
Hik which quantifies the cost for individual i to agree 
with individual k as 12311 



Hik{si,Sk) — —JikSiSk, 



(3) 



where Jik represents the strength of interaction between 
agents i and k. When i and k agree {siSk = 1) we have a 
cost Hik — — Jik, while when they disagree (siSk — — 1) 
we have Hik = Jik- Thus, the interaction works in such 
a way that, when Jik > 0, then i and k tend to imitate 
each others assuming the same trait and vice versa when 
Jik < 0. The magnitude of Jik gives how important it is 
for i to agree or disagree with k [20| . 

For the whole population we have the total cost func- 
tion 



if(s,J) = ^i7,fe = -^J, 



(4) 



(2) 



where the sum is extended over all the couples of nearest- 
neighbour agents denoted as k i. 

The cost function of Eq. |4]is just the well-known Ising 
Hamiltonian (see e.g. [^]) which can be treated by sta- 
tistical tools. As it is well known, the cost function 
-ff(s,J) does not lead to any natural dynamics and it 
is a very interesting matter of investigation to define 
a proper dynamics able to make the system evolve to- 
wards an "equilibrium" state. This can be achieved in 
several ways: Apart from exact analytical approaches, 
available only for special structure topologies (one di- 
mensional and two dimensional lattices) and mean-field 
solutions, a number of approximate techniques have been 
developed, including series expansions, field theoretical 
methods and computational methods. 

Here we adopt Monte Carlo (MC) numerical tech- 
niques in order to simulate the evolution of the system 
from a given initial configuration Sq to the stationary 
state, which in general does depend on the evolutionary 
dynamics D we choose and on the parameters Jik- 

Notice that J can be chosen to be directed and group 
dependent and this may account, in the examples dis- 
cussed in the introduction, for different influences and 
fares inter and intra different groups. For instance, if, 
say, the A company applies very high costs for phone 
calls between different users, we expect the pertaining in- 
teraction strengths to be very large. On the other hand, 
if for a B user costs for phone calls towards A users are 
only slightly more expensive then the relevant interaction 
strengths are small. 
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III. DIFFUSIVE STRATEGIC DYNAMICS 

In order to simulate the evolution of the population 
described by the cost function in Eq. [3J several different 
algorithms have been introduced. Among them a well- 
established one is the so-called single-flip algorithm which 
makes the system evolve by means of successive opinion- 
flips, where we call "flip" on the node j the transforma- 
tion Sj — Sj. 

More precisely, the algorithm is made up of two parts: 
first we need a rule according to which select an agent 
to be updated, then we need a probability distribution 
which states how likely the opinion-flip is. 

As for the latter, we adopt the well-known Glauber 
probability: Given a configuration s, the probability for 
the opinion-flip on the j-th node is given by 

p(s,j,J) = A . , (5) 



1 



,Ai/(s,j,J) 



where AiJ(s, j, J) is the variation in the cost function due 



to the flip Sj — ?> 



Notice that, for single- flip dynamics. 



the cost variation AH consequent to an opinion-flip only 
depends on the opinion of a few agents, viz. the j-th 
one undergoing the flipping process and its aj nearest- 
neighbours. This can be shown by spelling out the cost 
function variation appearing in Eq. Q: 



Ai/(s,j,J) =2s, V J,js, 



3 -Jij^i- 



(6) 



Interestingly, as can be derived from Eq. [Sj each 
opinion-flip is the result of a stochastic process featur- 
ing a competition between an energetic and an entropic 
term: the lower the cost of the opinion-flip and the more 
likely its occurrence. The external parameter J tunes 
the probability for an energetically unfavourable event 
to happen: For very low values of J any event is equally 
likely to happen independently of the magnetic configu- 
ration, conversely, for high values of J, when the agent 
J is surrounded by agents sharing the same opinion, the 
flipping of Sj gets a rare event. 

As already recalled, the opinion-flip probabilities just 
described can determine a dynamics only after a pre- 
scription for updating the system has been introduced. 
In other words, we first need a selection rule according 
to which extract agents, then the opinion of the selected 
agent will be possibly updated according to p{s,j,J). 
There exist several different choices for the first proce- 
dure, ranging from purely random to deterministic. 

Now, the most popular algorithms select nodes to be 
updated according to a sequential order which, though 
computationally efficient, appears rather artificial in a 
social network. Indeed, unless no predetermined strate- 
gies are at work, the random updating {D = TZ) seems to 
be the most plausible. In this case the probability that 
the current configuration s changes into due to the fiip 



-Si, reads 



I 



The dynamics generated by has been intensively 
studied in the past (see e.g. [2ai) and it has been shown 
to lead the system to the usual equilibrium (canonical) 
distribution, derived from the cost function H . 

However, in a social context, we notice that an opinion- 
fiip can occur as a result of a direct interaction (phone 
call, mail exchange, etc.) between two neighbours and if 
agent i has just undergone an opinion-fiip he will, in turn, 
pr omp t one out of his neighbours to change opinion 
[9l-ll2l|. This kind of picture can not be described by a 
random updating rule. Moreover, in many situations, the 
arbitrary agent i aims to be surrounded by neighbours j 
sharing his own opinion (being this a cultural trait or a 
phone subscription), i.e. SiSj — 1, because this translates 
in some advantage for agent i. 

Here we want to explore a relaxation dynamics, D — S, 
able to take into account these aspects, namely a realis- 
tic selection rule and a proper strategy. With respect to 
traditional dynamics, S displays two important features: 
i. the selection rule exhibits a diffusive character: The 
sequence of sites selected for the updating can be thought 
of as the path of a random walk moving on the graph rep- 
resenting the social network; ii. the diffusion is biased: 
The ai neighbours are not equally likely to be chosen 
but, amongst the at neighbours, the most likely to be se- 
lected is also the most likely to undergo an opinion-fiip, 
namely the one which minimizes AH{s, j, J). This corre- 
sponds to a local strategic choice of agent i, as the chance 
to obtain an opinion fiip is high, though the stochastic 
character is preserved. 

Let us now formalize how S works. Our MC simu- 
lations are made up of successive steps, each of them 
follows as: 

- Being i the newest updated agent (at the very first 
step i is extracted randomly from the whole set of 
agents), we consider the corresponding set of nearest- 
neighbours defined as A/i = {ii, 12, ia;}; we possibly 
consider also the subset Mi C Mi whose elements are 
nearest-neighbours of i not sharing the same opinion: 
j € M j € Mi A SiSj = —I. Now, for any j g Mi 
we compute the cost function variation A_ff(s,j,J), see 
Eq. [HI which would result if the flip Sj —Sj occurred; 
notice that AiJ(s,j, J) involves not only the nearest- 
neighbours of i, but also its next-nearest-neighbours. 

We calculate the probability of opinion- 
flip for all the nodes in A/i, hence obtaining 
p(s,ii, J),p(s,Z2, J), ...,p(s,ia,, J), where p(s,j,J), see 
Eq. O is the probability that the current configuration s 
changes due to a flip on the j-th site. 

- We calculate the probability V'^ {s;i, j; J) that, 
among all possible ai opinion-flips considered, just the 
j-th one is realized. This probability is defined as: 



_ p(s,j,J) 



fceA/i 



(8) 



7'^(s,j;J) = -p(s,j,J) 



(7) 



namely it follows by properly normalizing the p(s, j, J). 
Notice that, according to Eq. |51 among the nodes in- 
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eluded in TVi, those which are more hkely to be flipped 
are also the more likely to be selected for the spin-flip. 

We can possibly restrict the choice just to the set A/i, 
hence defining ^'^(s; J) = p(s, j, J)/ J^keSf. ^(s, k, J). 
The probabilities V'^ and V'^ are non trivially different 
from each other, due to the fact that they depend not 
only on the magnetic neighbourhood of the current site 
i, but also on the next-neighbourhood. As a consequence, 
in general, the next flip will not necessarily occur in A/i, 
not even for large coupling strength. Nonetheless, as 
we have verified numerically, the most striking result in- 
duced by the strategic dynamics, i.e. a shift of the critical 
region towards lower values of the interaction parameter 
(vedi infra) is preserved if we adopt instead of . 

- According to the normalized probability (see 
Eq. [8]), we extract randomly the node j G A/i and re- 
alize the opinion flip Sj — > ^Sj- 

- We set j = i and we iterate the procedure. 

Some remarks are in order now. The detailed balance 
condition usually implemented in standard dynamics and 
leading to a standard Boltzmann distribution, is explic- 
itly violated (see [1^ \2J} for more details). This does 
not contradict our dynamic intent; the evolution is not 
meant to recover a canonical Bolzmann equilibrium, but 
rather to model a realistic dynamics making the system 
evolve. By comparing Eq. [7] and Eq. |S1 we notice that 
the latter has an additional site dependence. As a re- 
sult, the analytical approach to the master equation is 
extremely difficult. Using the expressions for the proba- 
bilities in Eq.m the evolution equation for the probability 
P{s,i,t;J) that, at time t, the system is in the configu- 
ration s and the last spin flip occurred at site i reads: 

P(s,i,t+l;J)-P(s,i,<;J) = (9) 
V^{s; k, i- 3)P{s\, k, t- J) - r^{s- i, k; J)P(s, i, t; J), 

where s'i represents the spin conflguration which differs 
from s only for the value of Si. Interestingly, the ad- 
ditional dependence on the sites can be interpreted in 
term of a biased random walker, located on site i and 
implementing the dynamics by interacting with the spin 
system [2^. Therefore the master equation can be seen 
to describe, at a given time i, the evolution of the prob- 
ability for the spin variables to be in the configuration s 
and for the random walker to be located at site i. This 
makes the analytic approach difficult, even in the sim- 
plest lattices. Some general results, mainly numerical, 
have been obtained previously and in different contexts; 
we now briefly review some basic facts, while for more de- 
tailed information we refer the reader to [2^, [23| . Firstly, 
the diffusive dynamics is able to eventually drive the sys- 
tem to a thermodynamically well-behaved steady-state, 
characterized by the choice of the interaction coupling 
J and by the underlying topology, and also to recover 
the expected phase transition. However, the steady-state 
reached differs from the expected canonical equilibrium 
in a non-trivial way [26| . A particular evidence of this 



is provided by the fact that the critical temperature (or, 
similarly the critical coupling) obtained with the diffusive 
dynamics is signiflcantly larger than the one expected 
and such an effect can not be accounted for by a simple 
rescaling of the temperature |26j. Moreover, such results 
appear to be robust, as they hold also on general finite 
dimensional structures and for spin-1 models ^2^. 

Finally, it is worth underlying that such a diffusive dy- 
namics is intrinsically meant for finite systems. First of 
all this is consistent with our intent to model the dynam- 
ics within a community which is, indeed, finite. More- 
over, on finite, connected graphs without traps, a ran- 
dom walk is always recurrent and visits each site. This 
ensures that if our Monte Carlo simulation is run long 
enough, each agent is selected for the spin-flip a number 
of times sufficient to obtain a series of decorrelated states 
over which perform statistical averages. Interestingly, the 
rate at which a given site is visited depends non-trivially 
on its local magnetic environment as well as by the tem- 
perature and by the network topology. In particular, as 
shown numerically in [28|, [2§| , sites more likely to be se- 
lected for the updating correspond to borders between 
clusters. 



IV. NUMERICS 

As mentioned before, the analysis of the diffusive dy- 
namics has been carried out mainly from the numerical 
point of view by means of extensive Monte Carlo sim- 
ulations [1^. Here we focus on the particular case of 
interaction parameters Jik independent of the particular 
couple of agents considered, i.e. Jik = J, for any i,k; 
this allows to highlight the role of the dynamics leading 
the system to a stationary state and to understand how 
it possibly affects such stationary state and the average 
trait . 

In the simulation, once the network has been deflned, 
we place a binary variable Si on each node i and allow 
it to interact with its nearest-neighbors. Once the ex- 
ternal parameter J is flxed, the system is driven by the 
single- flip dynamics and it eventually relaxes to a station- 
ary state characterized by well-deflned properties. More 
precisely, for an ER random graph characterized by pa- 
rameters a and J, after a suitable time lapse to and for 
sufficiently large systems, measurements of a (specific) 
physical observable x{s,a,J) fluctuate around an aver- 
age value only depending on the external parameters J 
and a. 

We also verifled that, for a system of a given flnite size 
N, the extent of such fluctuations scales as N~2 (see 
also [l^ [13), as indicated by standard statistical me- 
chanics for a system in equilibrium. The estimate of the 
a given observable (x) is therefore obtained as an average 
over a suitable number of (uncorrelated) measurements 
performed when the system is reasonably close to the 
equilibrium regime. The estimate is further improved 
by averaging over different realizations of the underlying 



6 



6000 



5000 



4000 



3000 



2000 



1000 




1 p 
0.9 
0.8 
0.7 
0.6 

To.s 

0.4 
0.3 
0.2 
0.1 



Random Updating 



Diffusive 
Dynamics 



0.5 



1.5 



• AT = 6000,p = 0.0015 
X Af = 9000, p = 0.0010 



2 

J{z) 



2.5 



3.5 



FIG. 1: Thermalization of a ER random graph made up of 
= 6000 agents and a = 30. Two dramatically different ini- 
tial configurations are considered and compared: an ordered 
configuration with s,; = 1 for any i (black) and a random con- 
figuration with Si = 1 (si = — 1) with probability 1/2 (red). 



random graph with fixed a. In summary, 



1 

— Y^x(s{tr,)) 

n=l 



, tn^to+nT (10) 



where s{t) denotes the configuration of the system at 
time step t and T is the decorrelation parameter (i.e. 
the time, in units of spin flips, needed to decorrelate a 
given magnetic arrangement from the initial state); the 
symbol E denotes the average over different realizations 
of the graph. 

In general, during a MC run in a given sample we find sta- 
tistical errors which are significantly smaller than those 
arising from the ensemble averaging (see also [30]). 

We stress once again that the final state obtained with 
the diffusive dynamics is stable, well-defined and, in par- 
ticular, it does not depend on the initial conditions., i.e. 
it has all the properties of an equilibrium state. This is of 
course well-established for standard dynamics and it was 
also verified for our diffusive dynamics. An example is 
shown in Fig. 1 where, for a ER random graph, with given 
(a, J), the specific value around which m(t„) eventually 
fluctuates does not depend on the choice of the initial 
configuration selected for the simulation. To this aim we 
plotted TO(s(i„)) obtained starting with a completely or- 
dered configuration (toq = 1) smd with a completely ran- 
dom one. Moreover, we verified that the measurement 
of the (specific) observables within a macroscopic sub- 
system yields the same results characterizing the whole 
system. Similar results hold for scale free and small- world 
networks. 

In the following we focus on systems of sufficiently large 
size so to discard effects related to small N. For the ER 



FIG. 2: Average opinion (m) for a population of A = 6000 
(•) and A = 9000 (x) agents on a ER random graph with 
p = 0.0015 and p — 0.0010, respectively; the average number 
of nearest- neighbours is therefore the same for both systems, 
a = 9. Results obtained with a heath-bath dynamics (TZ) and 
with the strategic dynamics (5) are compared: for the latter 
a smaller critical parameter Jc is found. 



random graph, for a wide range of interaction constants 
J and average coordinations a, we measure the average 
magnetization {m)a.j (hereafter (m)) and the suscepti- 
bility X, calculated as 



xia,J) = JN [{m^)~{m)' 



(11) 



This quantity measures, at equilibrium, the reactivity of 
the system to a small external perturbation. Moreover, 
we compare results obtained for our dynamics with those 
obtained through a well-established algorithm, i.e. the 
Glauber heath-bath with random updating {D = TV), 
which is known to lead the system to a canonical steady 
state. 

In Fig. 2 we show results pertaining to ER random 
graphs of different sizes (A^ = 6000, A^ = 9000), but 
keeping the average coordination number fixed (a = 9 
corresponding to p = 0.0015 and p = 0.0010, respec- 
tively). Their profiles display the typical behaviour ex- 
pected for the Ising model on a random graph [Tsl] and, 
consistently with the theory, highlight a phase transition 
at a well defined value Jf (a)- Otherwise stated, there 
exists a critical value of the parameter J below which 
the system is spontaneously ordered. 

Note however that J^(a) is appreciably smaller than 
the critical value Jc{a) ~ 1/a expected for the canonical 
Ising model defined on the ER random graph [15J . 

An analogous behavior is also observed on scale-free 
and small-world networks; data for the average opinion 
(to) on these structures is shown in Fig. 3. More precisely, 
for the scale-free network considered here, the degree dis- 
tributions for nodes follows the power law pa with 
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FIG. 3: Average opinion (m) for a population of = 1000 
agents on a scale-free (left panel) and small-world (right 
panel) graph. Results obtained with a standard dynamics 
(TV) and the strategic dynamics (S) are compared: for the 
latter a smaller critical parameter Jc is found. 




10-^ 10-^ 



J 

FIG. 4: Average opinion (m) for a population of A'^ = 9000 
agents with a = 11 (•), a = 30 (x) and a = 45 (open circles) 
on a ER random graph. Results obtained with a standard 
dynamics (TZ) and the strategic dynamics (<S) are compared: 
for the latter a smaller critical parameter Jc is found. 



7 = 2.3; as for the small-world network, this has been 
built according to the prescriptions explained above and 
in such a way that the average degree per node equals 4. 

Remarkably, similar diffusive dynamics have been 
shown to lead an analogous decrease of the critical in- 
teraction parameter on re gula r structures [l^ and also 
for the spin-1 Ising model |27h29| . In these cases it was 
proved that a simple rescaling of the interaction constant 
J can not account for the differences between results pro- 



duced by the diffusive dynamics and a heath-bath dy- 
namics. This feature constitutes a first signature of the 
fact that the equilibria generated by the diffusive dynam- 
ics are not governed by the Boltzmann distribution. 

As mentioned above, on an ER random graph the criti- 
cal value Jc depends on the system size and on the proba- 
bility p, through their product a, i.e. Jc ~ 1/a. In order 
to check if a similar behaviour also holds for the dynam- 
ics S, we now fix the size of the system and make a vary; 
results for N = 9000 with a = 11,30,45 are reported 
in Fig. 4. Indeed, also for jf , we evidence a monotonic 
increase with a, however, in order to establish the actual 
relation between J and a, further extensive simulations 
are necessary. 

Similar to what happens with the usual dynamics, the 
relaxation time needed to drive the system sufficiently 
close to the equilibrium situation is found to depend on 
the parameter J. More precisely, we experience the so 
called critical slowing down: the closer J to its critical 
value, the longer the relaxation time. 

We now turn to the susceptibility defined in Eq. [TT] 
results for the ER random graph are shown in Fig. 5. 
In the thermodynamic limit, at the critical point Jc, the 
susceptibility diverges, while for finite sizes the suscep- 
tibility is expected to display a peak at Jc- This kind 
of behaviour is found also when the diffusive dynamics 
is applied and x just peaks at jf . An important point 
is that the shape of the curve is not modified, indicating 
that the reaction of the system to an external perturba- 
tion is conserved, with respect to the usual equilibrium, 
in the vicinity of the critical point. The social system 
is therefore expected to behave in the same way. Hence, 
we have further evidence of the fact that the diffusive 
strategic dynamics recovers the phase transition typical 
of the Ising model, though providing a lower value for 
the critical interaction parameter. 



V. CONCLUSIONS AND PERSPECTIVES 

In this paper, we have introduced a dynamics for so- 
cial systems displaying diffusive and strategic character. 
This dynamics has been shown to relax the system to 
thermodynamically well-behaved steady states; in par- 
ticular, after a suitable time, the values of the global 
observables of the system display averages independent 
of the initial conditions (which can, at least, affect the 
orientation of the asymptotic arrangement). The mag- 
netization, representing the average trait reached by the 
system as a function of the interaction parameter J, fea- 
tures a transition at a value of J which is strictly lower 
than the one obtained with a non-strategic random choice 
for the opinion fiips, on the same network. In particular, 
on the ER random graph, also the shape of the suscep- 
tibility near the transition point is conserved, indicating 
that the reaction of the social system to a small external 
perturbation in the stationary state is not modified. 

This picture indicates that with a strategic (local) 
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FIG. 5: Average susceptibility x for ^ population of A'^ = 9000 
agents with a — 30 (x) and a — 45 (o) on a ER random 
graph obtained with the strategic dynamics. Notice that the 
function peaks at a value smaller than the critical Jc expected 
for standard dynamics (dotted lines). 

choice in opinion flips a full-consensus configuration is 
obtained for lower values of the interaction parameter J, 
namely it is "easier" to obtain a community with an ori- 
ented opinion. Differently stated, in a society where a 
given value of the interaction J is present, the number 
of people with an oriented opinion is higher if the equi- 



librium is reached by a strategic opinion flip. The shift 
effect on the critical parameter is a general feature of 
the strategic dynamics, and our results have been shown 
to hold on very different topologies, such as ER random 
graph, scale-free random networks with hubs and small- 
world networks 16]. Transition points in social systems 
lead to extreme sensitivity of the global parameter to in- 
teractions, so its position and its shift represent a key 
information in the understanding of the stationary state. 
It can be interesting to enhance them or to avoid them, 
depending on the meaning of the global parameter. 

Finally, our dynamics allows a number of generaliza- 
tions concerning, for instance, the number of possible cul- 
tural traits or competitive companies admitted (namely 
the magnitude of the spin variable) [13] , the possible pres- 
ence of an external magnetic field (representing the effect 
of external biases such as advertisement) or a more com- 
plicated set of constant parameters Jij. In particular, J 
could be a directed, block matrix and this would account 
for different fares inter and intra distinct groups; in this 
case it would be very interesting to understand the con- 
ditions, in terms of J elements, for the realization of an 
oriented (i.e. magnetized) system. 
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